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ABSTRACT 



The existence of big bang relic neutrinos — exact analogues of the big bang relic 
photons comprising the cosmic microwave background radiation — is a basic pre- 
diction of standard cosmology. At present, the observational evidence for their 
existence rests entirely on cosmological measurements, such as the light elemen- 
tal abundances, anisotropies in the cosmic microwave background, and the large- 
^ . scale matter power spectrum. In this review, we concentrate on the prospects 

of more direct, weak interaction based relic neutrino detection techniques which 
. are sensitive to the cosmic neutrino background near the present epoch and in 

I/"") ' our local neighborhood in the universe. In this connection, we emphasize the 

. necessity to take into account the gravitational clustering of the relic neutrinos 

' in the cold dark matter halos of our Milky Way and nearby galaxy clusters. 

1. Introduction 

Over the last few years, significant advances were made in observational cosmolo- 
gy^. The position of the first Doppler peak in recent measurements of the cosmic mi- 
crowave background (CMB) radiation strongly suggests that the universe is spatially 
flat. Observations of Type la supernovae (SN) and the large scale structure (LSS) 
favor a universe with present energy fractions of 70 % of dark energy — accounting 
for the observed accelerating expansion of the universe — and about 25 % of cold 
dark matter — the corresponding particles being non-relativistic already at the time 
of recombination. Successful big-bang nucleosynthesis (BBN) of the light elements 
requires that about 5 % of the energy content of the universe is in the form of ordinary 
baryonic matter. All in all, we have now a pretty good knowledge of the cosmic recipe 
(cf. Table Q). A lot of further observational and theoretical effort will go into this field 
to further substantiate and explain these cosmological findings, and, in particular, to 
clarify the nature of dark energy and dark matter. 

Along with the CMB, standard big bang theory predicts the existence of a cosmic 
neutrino background (CuB). Presently, the evidence for the existence of the relic neu- 
trinos rests on the aforementioned cosmological measurements, i.e. the ligh t ele mental 
abundances, CMB anisotropies, and the large-scale matter power spectrunPEJ. Note, 

a Talk at XI International Workshop on "Neutrino Telescopes", Feb 22-25, 2005, Venice, Italy. 



Tabic 1: The cosmic recipe. 
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however, that all these measurements probe only the presence of the relic neutrinos 
at early stages in the cosmological evolution, and this often in a rather indirect way. 
It is therefore natural to ask: what are the prospects of a more direct, weak in- 
teraction based relic neutrino detection, sensitive in particular to the CzaB in the 
present epochs'? After all, among the known elementary particles, neutrinos are one 
of the most abundant particles in the present universe, falling second only to the relic 
photons (cf. Table HJ). 



2. How Many? How Fast? 

In order to design a direct, weak interaction based detection experiment, a precise 
knowledge of the phase space distribution of the relic neutrinos is indispensable. In 
this section, we will review recent determinations of this distribution"^*. 

The big bang relic neutrinos originate from the decoupling of the weak interac- 
tions when the universe was about one second old and the primordial plasma had 
a temperature of about one MeV, much larger than the possible neutrino masses. 
Neglecting late-time, small-scale gravitational clustering, to which we come later, 
and in the absence of appreciable lepton asymmetries, 13 their phase space distribu- 
tion is therefore predicted to be given by the homogenous and isotropic relativistic 
Fermi-Dirac distribution, fo(p) = 1/(1 + expfjo/T^o)), where p is the modulus of the 
comoving three-momentum p and T v $ = (4/11) 1 / 3 T 7i o = 1.95 K is today's neutrino 
temperature. Correspondingly, the gross properties of the Cz/B are tightly related to 

b Large neutrino m ixing infe rred from oscillation experiments ensures the validity of the neglection 
of chemical potential J7 I S I 9 1 1 ( J . 



0.1 



> 

<U 10" 2 



10- 3 r 



10- 



~~ 1 — 1 1 1 r 1 1 r J 1 — 1 1 1 [ 1 [ 1 J 1 — 1 1 1 t 1 r 1 1 

Normal hierarchy 
m. 




1 1 1 1 



10- 4 10- 3 10" 2 0.1 



0.1 



> 

V lO" 2 



10" 3 r 



10- 



1 1 1 1 1 1 1 1 — 1 — 1 1 1 1 1 1 1 1 — 1 — 1 1 1 1 in 
Inverted hierarchy 

m„ , m„ 




I I 1 I 



10- 4 10- 3 10- 2 0.1 



rn [eV] 



m„. [eV] 



Figure 1: Present knowledge about the neutrino mass spectrum^*!) , as a function of the mass of 
the lightest neutrino, m vi . In the normal hierarchical spectrum (left), the smaller mass difference 
inferred from solar neutrino oscillations separates the lightest neutrino from the next-to-lightest, 
whereas the larger mass difference inferred from atmospheric neutrino oscillations separates the 
latter from the heaviest neutrino. In the inverted hierarchical spectrum (right) this is reversed. 



the properties of the well-measured CMB and are therefore to be considered as rather 
firm predictions. Their present number density, 

3 

n Vi o = ^0 = — n-yo = 56 cm' 3 , (1) 
„ ' 22 

OB CMB 

when summed over all neutrino types i — 1,2, 3, is large and comparable to the one 
of the CMB, J2A^u t o + ^0) = (9/11) n 70 . Their present average three-momentum, 
on the other hand, is very small, 

P, i0 = Pp l0 = 3(^) 1/3 T 70 = 5x 10- 4 eV. (2) 

OB CMB 

Correspondingly, at least two of the relic neutrino mass eigenstates are non-relativistic 
today (m Ui 3> p^o), independently of whether neutrinos masses have a normal hier- 
archical or inverted hierarchical pattern (cf. Fig. HJ. These neutrinos are subject to 
gravitational clustering into gravitational potential wells due to existing CDM and 
baryonic structures, possibly causing the local neutrino number density to be en- 
hanced relative to the standard value (0) and the momentum distribution to deviate 
from the relativistic Fermi-Dirac distribution. 

This question can be studied quantitatively as followsPEI. First of all, in the 
context of a flat ACDM model, the neutrino contribution to the total dark matter 




Figure 2: CDM mass density profiles inferred from iV-body simulations of ACDM mo del . 

density is always a small perturbation (cf. Table Correspondingly, the CDM mass 
density p m dominates in the gravitational potential 0. Secondly, the neutrinos will 
have negligible gravitational interactions with each other. Therefore, one can make 
use of the CDM halo profiles from high-quality, pure ACDM simulations (cf. Fig. |5J) 
and study the evolution of the neutrino phase space distribution /^(x, p, r) in the 
corresponding gravitational potential wells. This distribution depends on the comov- 
ing distance x = r/a(t), its conjugate momentum p = am„. x, and the conformal 
time dr = dt/a(t), a being the cosmological scale factor. Technically speaking, one 
has then to solve the Vlasov, or collisionless Boltzmann equation, 

_ df Vi . df Ui df Vi . . 

Dr or ax. w ' ap 

p 

with the Poisson equation 

V 2 = 4vrGa 2 (p rn (x, r) - pjr) ) (4) 

V v ' 

<5 m (x,r)p m (r) 

relating the gravitational potential cf) to the CDM density fluctuation 5 m with respect 
to the physical mean p m . For a given CDM halo profile, e.g. the one advocated by 
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Figure 3: Neutrino number density profiles, normalized to their cosmological mea rP, calculated 
with ./V-one-body simulations (solid) and with the linear approximation (dotted), respectively. 
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where the parameters r s and p s are basically determined by the halo's virial mass 
M V i r , the Vlasov equation (JHJ) may be solved numerically by tracking the trajectories 
of neutrinos in phase space, starting from initial conditions corresponding to the ho- 
mogeneous and isotropic Fermi-Dirac distribution (iV-one-body simulations) - '. The 
initial redshift can be taken as z = 3, since, at higher redshifts, a sub-eV neutrino 
has too much thermal velocity to cluster efficiently 

A comprensive and exhaustive comparative study has revealed the following re- 
sults (cf. Fig. EI). First of all, a flattening of the neutrino profiles at small radii is 
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Figure 4: Neutrino density profiles for the Milky Waj^ , obtained via iV-one-body simulations from 
the MWnow (top curve in each plot) and the NFWhalo run (bottom curve in each plot). 

observed, in distinction to the CDM profiles (cf. Fig- EJ) ) - This can be understood in 
terms of neutrino free-streaming. Secondly, the clustering is considerably improved 
towards larger CDM halo virial masses and/or larger neutrino masses. In the inner 
part (<100 kpc) of a massive galaxy cluster like e.g. the nearby (~ 15 Mpc) Virgo 
cluster (M vir ~ 10 15 M Q ), the relic neutrino density n u can be larger than its cosmo- 
logical mean n u by a factor of ~ 10 4 1000, if we take into account the full range of 
possible neutrino masses for the heaviest neutrinos, m v = 0.04 4 0.6 eV (cf. Fig. 
Much more moderate clustering, n v /n v ~ 1 4- 20, is obtained for ordinary galaxies 
like the Milky Way (M vir ~ 1O 12 M ) in their central region (<10 kpc). We note 

in passing that a previously^ exploited semi-analytical linear method for solving the 
Vlasov equation (jSJ), which consists of replacing df/dp by df /dp, systematically 
underestimates the neutrino overdensities over nearly the whole range of halo and 
neutrino masses considered here (cf. Fig. EJ). 
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Figure 5: Momentum distribution of the relic neutrinos in the local neig hbourhood of the EartrP, 
obtained from the MWnow (solid) and NFWhalo (dashed) run, approaching for large momenta the 
relativistic Fermi-Dirac distribution (dotted). 



For an accurate determination of the relic neutrino number density and momen- 
tum distribution in the Earth's local neighbourhood in the Milky Way, at a distance 
r ffi ~ 8 kpc from the galactic center, one needs, in principle, to know the gravitational 
potential over the history of the Milky Way, i.e. the complete assembly history. In the 
absence of that information, one may consider two extreme cases, with the true be- 
havior som ehow l ying in-between: (i) working with the present day Milky Way mass 
distribution 1 15 1 1 ^ (MWnow), assuming it to be static (in physical coordinates), and 
(ii) exploiting the NFW halo (NFWhalo) that would have been there, had baryon 
compression — which is thought to lead to the formation of the galactic bulge and 
disk — not taken place. The possible ranges of overdensities, ~ 1 -v- 20, in the Milky 
Way are illustrated in Fig. HJ 

The final momentum distribution at r e is found to be almost isotropic, with mean 
radial velocity (v r ) ~ and second moments that satisfy approximately the relation 



Table 2: Relic neutrino properties^ as relevant for flux detection. 
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2(v 2 ) w {vj). The coarse-grained phase space densities f(r 9 ,p) in Fig. El are flat at 
low momenta, with a common value of nearly 1/2, have a turning point at about 
the escape momenta p esc = m u v esc = m„^2|0(r e )|, and quickly approach the Fermi- 
Dirac distribution for larger momenta. Note, that the results displayed in Fig. El 
n °t o nly sat isfy but, up to p esc , nearly completely saturate the general phase space 
bound^^ / < max(/o) — 1/2- The corresponding semi-degenerate state can only 
be made denser by filling in states above jt w. In order to attain even higher densities, 
one must appeal to non-standard theories^*. 



3. How to Detect? 



The gravitational infall of the relic neutrinos into CDM halos, discussed in the last 
section, might be significant for its influence on the non-lin ear (i.e. not so large-scale, 
wave number k>l Mpc -1 ) matter power spectrurrJ^^ , which will be determined 
in the upcoming weak gravitational lensing campaigns. In contrast to the evidences 
for the relic neutrinos from BBN, CMB, and the linear (i.e. large-scale, k<l Mpc -1 ) 
matter power spectrum mentioned in the Introduction, this gravitational inference 
will test the presence of relic neutrinos at low redshift, i.e. near to the present epoch. 
In this section, we will concentrate on other detection techniques, which are also 
sensitive to the present Cz/B, which, however, are based instead on weak interaction 
scattering processes involving the relic neutrinos either as a beam or as a targe P. 



3.1. Flux detection 



The Earth is moving through the almost isotropic (cf. last section) relic neu- 
trino background. In this subsection, we consider coherent elastic scattering of the 
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corresponding relic neutrino flux off target matter in a terrestrial detector. In this 
connection, it is an important observation that the low average momentum of relic 
neutrinos corresponds to a de Broglie wavelength of macroscopic dimension, 

* = l/(p> = 0.12 cm/(p/T ufi ). (6) 

Correspondingly, one may envisage^EHJ scattering processes in which many target 
atoms of atomic mass A act coherently over a macroscopic volume X 3 , yielding an 
enhancement of the elastic scattering rate by the huge factor 
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p t k 3 ~ 6 x 10 




(7) 



where Na is Avogadro's number and p t is the target mass density, compared to case 
where neutrinos are elastically scattered coherently only on the individual nuclei of 
the target. A terrestrial t arget of li near size r t < X will therefore experience a neutrino 
wind induced acceleratior P^EB ED 
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where a u ^ ~ G 2 F ml/ii is the elastic neutrino-nucleon cross section, and v re \ — (|v — 
v®|) is the mean velocity of the relic neutrinos in the rest system of the detector. 
Here, t>® ~ 7.7 x 10~ 4 c denotes the velocity of the Earth through the Milky Way. For 
Majorana neutrinos, the acceleration is further suppressed, in comparison with (jSJ), 
by a factor of (f re i/c) 2 — 10~ 6 for an unpolarized and v m \/c ~ 10 -3 for a polarized 
target, respectively. 



Table 3: Planned and projected accelerator beams and their interaction rates with the relic 
neutrinoJ^l . 
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What are the prospects to measure such small accelerations? Presently, con- 
ventional Cavendish-type torsion balances routinely reach 10~ 13 cm s~ 2 . Possible 
improvements with currently available technology to a sensitivity of > 10~ 23 cm s -2 
have been proposecP^^ (cf. Fig. |)J) C . However, even such an improved sensitivity 
is still off the prediction (JHJ) by at least three orders of magnitude, as an inspection 
of the currently allowed range of local relic neutrino overdensities in Table E] reveals. 
Therefore, an observation of this effect will not be possible within the upcoming 
decade. But it can still be envisaged in the not- so- distant future, say, within thirty to 
forty years. Note, in this context, that the acceleration can be improved still by a 
considerable amount by using foam-likd^ or laminated^ materials. In this way one 
may exploit a target size much larger than X, while still avoiding destructive inter- 
ference. Alternatively, grains of siz e ~ X could be randomly embedded (with spacing 
~ X) in a low density host materiaP32HJ . Advances in nanotechnology may be very 
welcome in this connection. For the case of Majorana neutrinos, flux detection via 
mechanical forces will remain a real challenge, however. 



3.2. Target detection 

The weak interaction cross sections are rapidly growing with energy, at least 
at center-of-mass energies below the W- and Z-resonances. In this subsection, we 
study the question whether the scattering of extremely energetic particles (accelerator 
beams or cosmic rays) off the relic neutrinos as a target has promising prospects for 
CvB detection. 

c Such improvements would also be very interesting in the connection of the search for deviations 
from Newton's law at small distances and possible significant improvements on the bounds of the 
size of extra space- like dimensional' . 



Figure 7: Resonant annihilation of an extremely energetic cosmic neutrino off a relic anti-neutrino 
into a Z-bosorP"' . 



3.2.1. Exploiting accelerator beams 

For a beam of particles %N, with energy E^, char ge Z e, len gth L, and current /, 
the interaction rate with the relic neutrinos is 



LI/(Ze) (9) 

* Tx 10- yr- (*L) (^) ^ 

y \nj \eVJ Z U0 TeV/ 



100 km/ V0.1 A 

In view of the currently allowed range of local relic neutrino overdensities, ~ 1 -j- 20, 
and the beam parameters of the next and next-to-next generation of accelerators — 
the Large Hadron Collider (LHC) and the Very Large Hadron Collider (VLHC), 
respectively — the expected rate (JSJ) is clearly too small to give rise to an observable 
effect in the foreseeable future (cf. Table 0J). Even with an Ultimate Large Hadron 
Collider (ULHC) designed to accelerate protons to energies above 10 7 TeV in a ring 
of ultimate circumference L ~ 4 x 10 4 km around the Earth, d it seems very difficult 
to establish the interactions with the relic neutrinos, although they occur at a rate of 
more than one event per year (cf. Table El): elastic scattering of the beam particles 
with the relic neutrinos — one of the contributions to the rate © — will be next to 
impossible to detect because of the small momentum transfers involved (~ 1 GeV at 



d Note that such an accelerator, in the collider mode, will probe the "intermediate" scale 
(Mew Mgvt) 1 ^ 2 ~ 10 10 GeV between the electroweak scale Mew ~ 1 TeV and the scale of grand 
unification Mqut ~ 10 17 GeV. The intermediate scale is exploited in many schemes of supersym- 
metry breaking and in seesaw mechanisms for neutrino masses. 




Figure 8: Signatures of resonant annihilation of extremely energetic cosmic neutrinos off relic neu- 
trinos: absorption dips (left) in the EECV flux^', for an isotropic source distribution out to redshift 
z max = 2, 5, 10 ( from upper to lower curves), and emission features (Z-bursts) (right) in the EEC ray 
(EECR) fkul^al flux in the form of photons (short-dashed) and protons (dashed-dotted) , eventually 
overcoming the ordinary extragalactic EECR flux (long dashed) , which suffers from the GZK cutoff. 



En ~ 10 7 TeV). A very promising alternative is to consider a heavy ion beam, and 
to exploit the contribution of the inverse beta decay reaction, 

A z N + v e ^ A z+l N + e- , (10) 

to the rate (jOJ) . This reaction changes the charge of the nucleus, causing it to follow an 
extraordinary tr aject ory and finally to exit the machine such that it becomes suscep- 
tible to detection JU | J A detection of this reaction would also clearly demonstrate 
that a neutrino was involved in the scattering. 



3.2.2. Exploiting cosmic rays 

In the foreseeable future, before the commissioning of the ULHC, target detection 
of the relic neutrinos has to rely on extremely energetic cosmic rays. Indeed, the 
resonant annihilation of extremely energetic cosmic neutrinos (EECz^s) off relic anti- 
neutrinos into Z-bosons (cf. Fig. EJ), occuring at the resonance energies 

ET = <_ „ 4 x 10 2i ( eV 

offers unique opportuniti es for relic neutrino detection. On the one hand, one may 
search for absorption dip d * ' * i \ i.L \ oo \ o\) \ j n ^ e EECz/ spectrum at the resonant 
en ergies (cf. Fig. [ H] (left)), on the other hand, one may look for emission featu- 
re j4U | 41 | 42 | 43 | 44 ^31 ^Z-bursts) as protons or photons with energies s pannin g a decade 
or more above the predicted Greisen-Zatsepin-Kuzmin (GZK) cutoflP^^ at Eqzk — 
4 x 10 19 eV (cf. Fig. |H1 (right)). This is the energy beyond which the CMB is absorbing 
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Figure 9: EEC, fluxeJ 1 M^NI^m 

current limits and projected sensitivities of dedicated 
experiments (left), as well as prospects within the next decade (2008 and 2013) and theoretical 
predictions (right). 
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Figure 10: Relic neutrino absorption dips in the EEO fhrjnj from hidden topological defect sources 
(left) and topological defect sources (right), with projected error bars in the year 2013. 



to nucleons due to resonant photopion production. Indeed, the association of Z- 
bursts with the mysterious post-GZK cosm ic rays ob served^" by the Akeno G iant 
Air Shower Array (AGASA) is a controversia 
Fig. El (right)). 

Presently planned neutrino detectors (Pierre Auger Observatory®!, IceCube^^, 
ANITA^l, EUSO^, OWlP3, and SalSA^S) operating in the energy regime above 
10 21 eV appear to be sensitive enough to lead us, within the next decade (cf. Fig. |0J), 
into an era of relic neutrino absorption spectroscopy ^ H^ l ^ t (cf. Fig. 1X0)1. provided 
that the EECu flux at the resonant energies is close to current observational bounds 
and the neutrino mass is sufficiently large, m v > 0.1 eV. In this context it is important 
to note, that absorption spectroscopy is predominantly sensitive to the relic neutrino 
background at early times, with the depths of the absorption dips determined largely 
by the higher number densities at large redshifts (z ^> 1) (cf. Fig. El (left)). Since 
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Figure 11: "Large scale" overdensities (i — v, CDM) in the local universe^, with the Milky Way 
at the origin. The black (solid) line corresponds to the local CDM distribution^^ , inferred from 
peculiar velocity measurements^' , smeared over the surface of a sphere with radius r. The dotted 
line is the neutrino overdensity for m v = 0.6 eV, short dash 0.3 eV, long dash 0.15 eV, and dot-dash 
0.04 eV. 

neutrinos do not cluster significantly until after z<2, clustering at recent times can 
only show up as secondary dips with such minimal widths in energjl^ that they do 
not seem likely to be resolved by planned observatories. 

On the other hand, emission spectroscopy is directly sensitive to the relic neutrino 
content of the local universe (z<0.01 ?"gzk<50 Mpc). However, since the neu- 
trino density contrasts approximately track those of the underlying CDM above the 
neutrino free-streaming scale, it is clear that there cannot be a substantial neutrino 
overdensity over the whole GZK volume (~ ^gzk)- Indeed, the estimated neutrino 
overdensity in our local GZK zone, with a ~ 5 Mpc smoothing, is always <2 (cf. 
Fig. 111)). Hence, the overall emission rate cannot be significantly enhanced by gravi- 
tational clustering. 

Nevertheless, it seems wort hwhile to contemplate about "relic neutrino tomogra- 
phy" of the local universe ^^* . Specifically, one may exploit the fact that there are 
several galaxy clusters (> 1O 14 M ), such as Virgo (distance ~ 15 Mpc) and Centaurus 
(~ 45 Mpc), within the GZK zone (cf. Figure IT2|) within which we expect significant 
neutrino clustering (cf. Figure EJ). One could then search for directional dependences 
in the emission events as a signature of EECu annihilating on relic anti-neutrinos 
(and vice versaf^. For example, the angular resolution of AGASA, ~ 2°, is already 
sufficient to resolve the internal structures of, say, the Virgo cluster (M vir ~ 1O 15 M ) 
which spans some 10° across the sky. Using the iV-one-body clustering results in 
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Figure 12: Various massive galaxy superclustera"^ in our "vicinity". The Milky Way is at the center 
of the coordinate system. 

Figure El the average neutrino overdensity along the line of sight towards and up to 
Virgo is estimated to be ~ 45 and ~ 5 for m v = 0.6 eV and 0.15 eV respectively, 
given an angular resolution of ~ 2°. The corresponding increases in the number of 
events coming from the direction of the Virgo cluster relative to the unclustered case, 
assuming an isotropic distribution of EECu sources, are given roughly by the same 
numbers, since protons originating from ~ 15 Mpc away arrive at Earth approxi- 
mately unattenuated. The numbers improve to ~ 55 and ~ 8 respectively with a 
finer ~ 1° angular resolution. In the most optimistic case of a EECu flux near to 
the current observational bound (cf. Fig. El (right)) and a neutrinos mass >0.1 eV, 
EUSO will not only find evidence for the absorption dips, but also for the enhanced 



emission from the direction of Virgo due to Z-burst; 



4. Conclusions 

At present, BBN, CMB, and the large-scale matter power spectrum provide the 
only observational evidence for the big bang relic neutrinos, at least in the early 
stages of the cosmological evolution. A more direct, weak interaction based detection 
of the Cz/B near the present epoch may proceed in the following chronological order 
by measuring 

i) absorption dips in EECV spectra and Z-bursts in EECR spectra; 

ii) macroscopic forces through coherent elastic scattering of relic neutrinos off tar- 
get material in Cavendish-type torsion balances; 

iii) interactions of extremely energetic particles from terrestrial accelerator beams 
with the relic neutrinos as a target. 

Unfortunately, an immediate and guaranteed direct detection does not appear to 
be feasible. Although the search for signatures of EECu annihilation off the relic 
neutrinos can start right now, its success entirely rests on the existence of an EECz/ 
flux at the resonance energies. The sensitivity of Cavendish-type torsion balances has 
still to be improved by at least thirteen orders of magnitude for a detectable signal, 
which postpones this detection techniques probably beyond the year of retirement 
of the author, >2025. Finally, an appreciable rate of beam particles with the relic 
neutrinos requires an Ultimate Large Hadron Collider around the Earth with a beam 
energy E^ eam > 10 7 TeV, which almost certainly will never be built. In the meantime, 
we can hope to see the late-time relic neutrinos through their gravitational effects 
in the not-so-large-scale, non-linear part of the matter power spectrum measured by 
weak gravitational lensing. 
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